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We study the real time formalism of non-equilibrium many-body theory, in a first 
quantised language. We argue that on quantising the relativistic scalar particle in 
spacetime with Minkowski signature, we should study both propagations e l<yV ~ m ^ and 
e -t(p -m )a on particle world line. The path integral needs regulation at the mass 
shell p 2 = m 2 . If we regulate the two propagations independently we get the Feynman 
propagator in the vacuum, and its complex conjugate. But if the regulation mixes the two 
propagations then we get the matrix propagator appropriate to perturbation theory in a 
particle flux. This formalism unifies the special cases of thermal fluxes in flat space and 
the fluxes 'created' by Cosmological expansion, and also gives covariance under change of 
particle definition in curved space. We comment briefly on the proposed application to 
closed strings, where we argue that coherent fields and 'exponential of quadratic' particle 
fluxes must both be used to define the background for perturbation theory. 
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1. Introduction. 

When we quantise a field in curved space, the notion of a particle becomes very 
curious. For example in an expanding Universe we have particle creation, which means 
that a spacetime which looks empty in terms of particle modes natural in the past, may 
look full of particles on using co-ordinates natural to the future. The Minkowski vacuum 
appears to have a particle flux for an accelerated observer. For black holes, imposing 
vacuum conditions at past null infinity gives Hawking radiation in the future, due to the 
time-dependent gravitational field of a collapsing object. (JIJ and references therein.) 

How 'real' are such particles, and more specifically, how do they affect the gravita- 
tional field? Beyond semi-classical approximations, such questions have traditionally been 
deferred to a time when a consistent theory is available with both matter and gravity 
quantised. 

Strings provide such a consistent theory, so we wonder if the above questions have 
been implicitly answered in computing string amplitudes. One deals with closed strings 
(which include the graviton in their spectrum) in a first quantised language. The particle 
analogue of this approach is summing over all particle trajectories between initial and final 
points (with branching trajectories giving possible interactions). If we just take a free 
scalar particle and compute the two-point function by such a first quantised path integral, 
which notion of particle are we using? 

The answer to this question is known ||. If |0 >j n is the vacuum based on positive 
frequency modes at t = — oo and |0 > ou t for the modes based at t = oo then 

J D[paths(x -> x')}e lS = out < Q\T[<j>{x')<j>{x)]\Q > in / out < 0|0 > in (1.1) 

(The notion of 'in' and 'out' vacuua here does not necessarily require flat spacetime at 
t -> ±00.) 

But we may not want such a hybrid 'in-out' expectation value for our two-point 
function. We may wish to use the notion of 'in' particles, which requires computing 
i n < 0\T[(p(x')4>(x)} |0 >i n . Or we may have an ensemble of 'in' particles to start with, 
which implies a density matrix based on 'in' states (e.g. p ~ e~^ n \n >i n j n < n\) with 
propagator Tr{pT[(f)(x')(j)(x)}}. How do we modify the first quantised path integral to 
obtain the requisite propagators, and carry out perturbation theory? 

To handle non-equilibrium situations with general density matrices one uses the 'real 
time' formalism to develop a perturbation theory. If p is specified in terms of the states 
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at t = — oo then we evolve the fields from t = — oo to t = oo and back to — oo where 
we insert p and take a trace. The propagator becomes a 2x2 matrix propagator, with 
a = 1, 2 labelling operators on the first and second parts of the above time path. With 
correspondingly generalised interaction vertices, one computes Feynman diagrams in the 
usual way to obtain the correlators in the many-body situation (|| and references therein). 

The goal of this paper is to study this real-time formalism from a first quantised 
viewpoint, using the simple example of a scalar field. We require a covariant language for 
the density matrix (and time path) rather than a Hamiltonian language based on spacelike 
slices, because it is such a covariant language that we can extend to strings. 

In brief, our approach and results are as follows. For flat space, we know that a density 
matrix of the form 'exponential of linear in the field' gives a coherent state. Shifting the 
field by its classical value removes such a part of p. The class of p of the form 'exponential 
of quadratic in the field' is also special; for correlators with such p the Wick decomposition 
holds Q. The matrix propagator of the real-time formalism encodes such p and the choice 
of time path for the perturbation theory. Departures from this 'exponential of quadratic' 
form of the density matrix gives 'correlation kernels', which are handled perturbatively as 
vertices analogous to the interaction vertices in the Lagrangian. 

For the first quantised language, we argue that a careful quantisation of the relativistic 
scalar particle requires considering both the propagation e l ( p ~ m -* and the propagation 
e -i(p -m ) on ^Yiq world line. We may collect these two possibilities into a 2x2 matrix 
form, getting a world line Hamiltonian diag[— (p 2 — m 2 ), (p 2 — m )]. Near the mass shell 
p 2 — m 2 = the path integral needs regulation. If we add — iel to the Hamiltonian (/ 
is the 2x2 identity matrix) then we get a diagonal matrix propagator, with (|1 . 1|) and its 
complex conjugate as the first and second diagonal entries. But if we regulate instead 
by adding —ieM, where M is a non-diagonal matrix, then we get the matrix propagator 
corresponding to an 'exponential of quadratic' particle flux. 

We compute M for a thermal distribution in Minkowski space, for the Niemi-Semenoff 
U choice of time path, and for the 'closed time path' developed for non-equilibrium theory 
by Keldysh ||, Schwinger 0, and others. These two M matrices are not the same, which 
reflects the fact that the matrix propagator depends on both p and the choice of time path. 
For a curved space example, we consider a 1 + 1 spacetime with 'sudden' expansion. We 
compute M to obtain the matrix propagator appropriate to perturbation theory with the 
'in' vacuum; i.e. for p = i n \0 >< 0\i n and time path beginning and ending at past infinity. 
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Thus the 'exponential of quadratic' form of the density matrix, basic to perturbation 
theory, is naturally obtained in the first quantised formalism. Thermal fluxes in flat space, 
and the flux given by the Bogoliubov transformation due to spacetime expansion, are 
special cases of this form. Further, this class of density matrices is closed under change 
of the basis of functions used to define particles, so our formalism is covariant under such 
transformations. 

We conclude with a discussion of the significance of our results for strings, which were 
the motivation for this study of the first quantised formalism. The usual first quantised 
path integrals for strings would give the analogue of ( |1 . 1| ) , which corresponds to specific 
boundary conditions at spacetime infinity. To handle phenomena involving particle fluxes, 
we propose extending the world sheet Hamiltonian to a 2x2 Hamiltonian as in the above 
particle case. Thus we would allow not only classical deformations of the background 
field (analogous to the 'exponential of linear' p in the particle case) but also particle flux 
backgrounds (corresponding to 'exponential of quadratic' density matrices). Studying (3- 
function equations || for this extended theory should give not equations between classical 
fields but equations relating fields and fluxes. The latter kind of equation, we believe, 
would be natural to a description of quantised matter plus gravity. 

fThe plan of this paper is as follows. Section 2 reviews finite temperature perturbation 
theory, and discusses the significance of 'exponential of quadratic' density matrices for the 
curved space theory. Section 3 translates the finite temperature results of flat space to 
first quantised language. Section 4 gives a curved space example. Section 5 is a summary 
and a discussion relating to strings. 

2. Propagators in the presence of a particle flux. 

2.1. Review of the real time formalism. 

Consider a scalar field in Minkowski spacetime (metric signatureH — . . . — ) 



Suppose this field is at a temperature T = (/3) 1 . To take into account this temperature 
we evolve the field theory in imaginary time from t = to t = —ifl, and identify these 
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two time slices. In the 'real time' approach of Niemi and Semenoff we use instead the 
following path C\ in complex t space to connect these two slices 



C\\ I : — oo — ► oo 

II: oo -> oo-i/3/2 

(2.2) 

III: oo-i/3/2 -> -oo-i/3/2 
IV : -oo - z/3/2 -> -oo - i(3 

The concept of time ordered correlation functions is now replaced by 'path ordered' cor- 
relation functions, with t running along the above path C\. One argues that the parts II, 
IV of C\ can be ignored. Field operators on part / are labelled with the subscript 1 while 
on part II are labelled with the subscript 2. Thus for the two point function of the scalar 
field there are four possible combinations of subscripts, and the corresponding propagators 
are collected into a matrix: 

n , , _ ( p ,_^ +ie +2nn(p)5(p 2 -m 2 ) 2ix[{n(p){n(p) + l)Y' 2 5(p 2 - m 2 )\ 

P 2 



2ir[n(p)(n(p) + l)} 1/2 5(p 2 - m 2 ) _ a _ m l_ ie + 2nn(p)S(p 2 - m 2 ' 



where 

n(p) = (e^ -l)" 1 (2.4) 

is the number density of particles for the free scalar field at temperature (/3) _1 . We can 
write 

D(p) = U(p)D°(p)U(p) (2.5) 



with 



and 



D°(p) = ( pA -^ +ie " ) (2.6) 



u(v) _ (V 1 + n (P) V n (P) \ (27 ) 

Perturbation in the coupling A is described by the following diagram rules. The 
external insertions are all of type 1. The vertices from the perturbation term must have 
either all legs of type 1 or all legs of type 2. Attach a factor (— iX) to the vertex in 
the former case, (zA) in the latter. The propagator is given by the matrix ( |2.3|) and can 
connect vertices of type 1 to vertices of type 2 at nonzero temperature. Summing all 
Feynman diagrams with these rules gives the correlation function at temperature 
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The power of the real time approach is that the density matrix need not be thermal, 
and the system need not be in equilibrium. Suppose the density matrix is specified at 
t = — oo in terms of the free field operators, so that the perturbation Hamiltonian will 
modify the distribution as time progresses. One uses the closed time path Ci ||, @: 

C2 = I ■ —00 — > 00 

(2.8) 

II : 00 — ► —00 

for time evolution, and takes a trace after inserting p. Perturbation theory may be devel- 
oped with the same rules as above, but using the matrix propagator appropriate to the 
contour Ci. In particular, if p at t = —00 is of thermal form with temperature then 
the matrix propagator is 

,2 /r ^2\ n„\(„(„\ i_ Q( ^ 



— f + 27in(p)S(p 2 — m 2 ) 27r[(n(p) + 0(—p )]5(p 2 — m z 

—i 

~~ 5 5 



" " 1 2n[n(p) + 9(p )]5(p 2 -m 2 ) - + 2nn(p)5(p 2 - m 2 ) ] 



Perturbation theory based on the contours Ci and C2 are not equivalent, even for zero 
temperature |J . Contour C2 corresponds to computing correlators with p = |0 >< 0| 
inserted at t = —00. Consider the scalar field in + 1 spacetime dimensions, for simplicity. 
Let the only perturbation be the time dependent term p f dtS(t — to) : 4> 2 (t) :, and compute 
to first order the two point correlator for field insertions at t\ < to < t^- 



< (t>(t 2 )(t>(ti) > = 2p[< O|0(* 2 )0(t o )|O >< O|0(t o )0(ti)|O > 

+ < 01^)^2)10 X O|0(t o )0(ti)|O >] 



(2.10) 



With the contour C±, for (3 — > 00, we just get the first term on the RHS. The second term 
on the RHS comes from correcting the state at t = 00 away from the vacuum |0 >, to be 
such that it evolves from |0 > after suffering the perturbation at to- I n short, the time 
path is significant because perturbation theory works in the interaction picture, while the 
physical situation is described in the Heisenberg picture. The second leg of the time path 
makes the state at t = 00 the same as the state at t = —00 in the Heisenberg picture, 
which implies that perturbation corrections must be made along this path segment in the 
interaction picture. 

Note that the corrections to D°(p) (eq. ( |2.6| )) that give ( p73| ) or ( pD| ) are all on shell. 
This is a manifestation of the fact that these corrections arise from a flux of real particles. 
The expression 5(p 2 — m 2 ) does not have a Euclidean counterpart. (In Euclidean space one 
would consider (p 2 + m 2 ) which is positive definite.) Thus temperature is a phenomenon 
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of Minkowski signature spacetime, where Green's functions are subject to the addition of 
solutions of the homogeneous field equation. 

To understand the origin of the on shell terms, let us consider the element D\\ in 
( |2.3|) . We can decompose the free scalar field in Minkowski space into Fourier modes. Each 
quantised mode is a harmonic oscillator with frequency u(p) = (p 2 + m 2 ) 1 / 2 . The two 
point function < T[<fi(x 2 )(J)(xi)} > reduces, for each Fourier mode, to < T[q(t 2 )q(ti)} > for 
each harmonic oscillator. Thus we focus on a single oscillator (we suppress its momentum 
label). Let t 2 > t\. Then 

< q{t2)q(ti) > = J2 e ~ 0(n+h) " < n\T[q(t 2 )q(ti)]\n > /£V0(»+*)" 

(2.11) 

= < n+1 > e -w(t a -tQ , < n > c iu,(t 2 -ti) 
2co 2u 

The first term on the RHS corresponds to the stimulated emission of a quantum at t\ with 
absorption at t 2 . The second term corresponds to the annihilation at time t\ of one of the 
existing quanta in the thermal bath, and the subsequent transport of a hole from t\ to t 2 , 
where another particle is emitted to replace the one absorbed from the bath. The Fourier 
transform in time of (|2.11[ ) is 



poo ■ 

/ dte~ wt < T[g(* 2 )g(*i)] > = — 3 + 2n < n > 5(u' 2 - u 2 ) (2.12) 

7-oo ' uj' 2 -uj 2 + ie 



which gives the matrix element Dn in ( |2.3| ). 

Thus the correction to < T[(f)(x 2 )(f)(xi)] >= Dn(x 2 ,xi) (and other elements of D) 
due to the particle flux is not an effect of interactions with the particles of the bath. This 
correction arises from the possible exchange of the propagating particle with identical real 
particles in the ambient flux. Using such a corrected propagator with the interaction 
vertices gives the interaction of the bath particles with the propagating particle. 

2.2. 'Exponential of quadratic' density matrices. 

Perturbation theory in the vacuum involves separating a free part which is described 
by propagators, and an interaction part which is described by vertices. In studying kinetic 
theory in Minkowski space, we also need to identify a free part and an interaction term. 
The free part involves specifying a density matrix p of a special form, and a choice of time 
path Q . These special p are of the form 'exponential of an expression with quadratic and 
linear terms in the field' Q. The linear term deccribes a coherent state, and gives a change 
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in the classical value of the field. Shifting the field by this classical value gets rid of this 
linear term, and we will always assume that this has been done. The quadratic part of the 
exponential implies a density matrix of the form 

p = TT e «<aJoi e -/3ioJoi e 7iOiOi (2.13) 
i 

where the product is over different frequency modes, and a^, at are the annihilation and cre- 
ation operators for these modes. We will call the form (|2.13|) an 'exponential of quadratic' 
density matrix. The special role of density matrices ( |2.13|) is due to the fact that Wick's 
theorem extends to correlators computed with such p [|J. Thus for operators Ai linear in 
the field, 



< A x ...A n > = 7 ^-Tt{ p A 1 . . . A n } = < A h A i2 > . . . < A in _ x A in > 



permutations 

(2.14) 

We sketch a proof of (|2.14| ) in the appendix. 

The time path specifies where p directly gives the particle flux. Thus a time path 
beginning and ending at t = — oo says that the density matrix (specified in the interaction 
picture) is the flux at t = — oo, which implies that the flux at later times would be modified 
by perturbative corrections. 

Perturbation theory may be developed for p of the form ( |2.13|) in the same manner as 
for the vacuum theory, with the replacement of propagators with matrix propagators as 
in the thermal examples of sec. 2.1. The matrix propagator encodes both the choice of p 
and the choice of time path. Deviations from the form ( |2.13| ) of p give rise to 'correlation 
kernels' which are three and higher point vertices arising from correlations encoded in the 
density matrix, rather than the interaction Hamiltonian Q. 

We take the point of view that it is incorrect to construct a theory of quantised 
matter and gravity without including 'exponential of quadratic' density matrices in the 
possible backgrounds about which the perturbation will be developed. Suppose we choose 
a time co-ordinate t and start with a distribution e~@ H , thermal for the Hamiltonian 
giving evolution in t. As the Universe expands, the distribution will not remain thermal, 
in general. Redshifting of wavelengths gives an obvious departure from thermal form if the 
field has a mass or is not conformally coupled. But even a massless conformally coupled 
field departs from thermal form if the time co-ordinate t is not appropriately chosen ||1Q| . 
However, the density matrix remains within the class ( |2.13| ) , if it starts in this class, even 



for a massive field. Indeed, we can describe such p in a covariant fashion by choosing a 
set of global solutions to the wave-equation, and attaching operators Oj, a\ to pairs of 
functions fi, f*. p gives the linear map (on this space of solutions) that must be made 
before identifying the two ends of the time path of the perturbation theory. 

The class ( |2.13p of p is closed under change of the basis of functions used to define 
creation and annihilation operators. For example an expanding Universe suffers particle 
creation, so that an initial vacuum state would be seen as filled with particles by an 
observer using 'out' frequency modes [0. The density matrix in terms of 'in' modes is 
|0 >in in < 0| which is obtained from (|2.13|) with a = 5 = 0, (3 — > oo. In terms of 'out' 



modes the Bogoliubov transformation gives the form p ~ e batat |0 > ou t out < 0|e 6 aa which 



is ( |2.13[ ) with a = 6, 7 = 6* , (3 — > oo. The fact that in each case we are describing the 
state at the past time boundary is encoded in the time path, which starts at this boundary, 
and returns to it. 

To summarise, the class of 'exponential of quadratic' p is natural for defining the 
propagator. Considering this class unifies the special cases of thermal fluxes in flat space 
and the fluxes 'created' in spacetime expansion, and also gives covariance under change of 
the basis functions in spacetime used to define particles. 



3. The first quantised formalism. 

3.1. The regulator matrix. 

The Feynman propagator for a scalar field in Minkowski space can be written as 

poo 

D F (p) = — = / d\e lX ^ 2 - m2+ ^ (3.1) 



p2 _ m 2 _|_ le J\=0 

fluTTD can be used to express Gf(p) in a first quantised language, with p 2 = — □. (See for 
example |Tl], ||12|| .) The Hamiltonian on the world line is — (p 2 — m 2 ), evolution takes place 
in a fictitious time for a duration A, and this length A of the world line is summed over all 
values from to oo. 

Let us write the matrix propagator ( [2.3| ) in a similar fashion 

D NS (p) = / d\e- fxH -^ XM (3.2) 



where 

H= f-(p 2 -m 2 ) \ / 1 + 2n(p) -2Vn(p)(n(p) + l)\ 

V (p 2 -m 2 ))' \-2 v /n(p)(n(p) + l) 1 + 2n(p) J 

(3.3) 

is given by (|2.4|) .) This matrix world line Hamiltonian has the following structure. 
If we forget the term multiplying e, then the a = 1 component of the vector state on the 
world line evolves as e lX ( p ~ m ) while the a = 2 component evolves as e~ lX ^ p ~ m \ To 
define the path integral we need the regulation from e at the mass shell p 2 — m 2 = 0. 
But the regulator matrix is not diagonal, for nonzero temperature. Thus transitions are 
allowed from state 1 to state 2. In the limit e — ► + , which we must finally take, these 
transitions occur only on the mass shell. (The absolute values of the entries in M are not 
significant, because e goes to 0, but the relative values are.) 

Similarily we can write the matrix propagator (|2.9|) in the form (|3.2| ) with H as in 
(p|) but 

_ / l + 2ra(p) -2y/n{p){n{p) + l) e -^°/ 2 

V -2y/n(p)(n(p) + l)e^°/ 2 1 + 2n(p) 

Again we have the a = 1, 2 states propagating on the world line with Hamiltonians =F(p 2 — 
m 2 ). Transitions between a = 1, 2 are again of order e, but are different from those in 

(ID- 

Looking at these examples the following picture emerges. To obtain the matrix propa- 
gator in a many-body situation (with 'exponential of quadratic' density matrices) we need 
to consider both evolutions e~ lXH and e lXH on the world line. A damping factor is needed 
to define the first quantised path integral, near the mass shell. But the regulator matrix 
M need not be diagonal, and it reflects both the 'exponential of quadratic particle flux 
and the choice of time path. (For example, in the limit (3 —>■ oo M in ( |3.3[ ) is diagonal, but 
in (ETij) it is not.) 

3.2. Quantising the relativistic particle. 

What is the origin of the two components a = 1, 2 of the state on the world line? We 
would like to offer the following heuristic 'derivation' as a more physical description of the 
matrix structure in (|3.2|). 

The geometric action for a scalar particle is 





where r is an arbitrary parametrisation of the world line. The canonical momenta 



satisfy the constraints 



P - 91 - mX ^ T (3 6) 



P^P^-m 2 = (3.7) 



We choose the range of the parameter r as [0, 1]. Following the approach in we impose 
the constraint at each r through a 5-function: 

5(p 2 (r)-m 2 ) = j- d\(r)e-^ p2 ^- m ^ (3.8) 

J — oo 

The path integral amplitude to propagate from Xi to Xf becomes 

f D[X]D[P]D[X] j C dr[P M (r)X^, T (T)-A/2(r)(P 2 (r)-m 2 )] /o q"\ 

voipiff] e ° (3 - 9) 

where iV is a normalisation constant, P fI X^ 1T = m(X fI 1T X M , T ) 1 / 2 is the original action 
( ft.5[ ) and we have divided by the volume of the symmetry group, which which is related 
to r-diffeomorphisms in the manner discussed below. (The 5-function constraint on the 
momenta and dividing by Vol[Diff] remove the two phase space co-ordinates redundant in 
the description of the particle path.) 

There are two ways to consider the symmetry of the action ( |3.9| ). The action is 
invariant under 

51 : SX^(t) = h{r)P^(r) 

5P M (r) = (3.10) 
8\{t) = h(r), T 

and 

S 2 : 5X^(t) = e(r)X(r)P^(r) 

SP^t) = (3.11) 
8X(r) = (e(r)A(r)), T 

The difference between Si and 52 is best seen by considering the finite transformations on 
A: 

Si : A'(r) = A(r) + ^p- (3.12) 

or 

5 2 : A'(r) = ^(r)A(r) (3.13) 
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Using <Si we can gauge fix any function A(r) to any other function Ai(r), provided A, Ai 
have the same value of 

1 

A(r)dr = A (3.14) 

o 

With ^2, A transforms as an einbein under the diffeomorphism r — > t'(t). Note that 
for regular e(r), A either changes sign for no r or for all r. We take Diff as the group 
of regular diffeomorphisms connected to the identity; then we have only the former case. 
These diffeomorphisms cannot gauge-fix A(r) to any preassigned function Ai(r). We again 
have the restriction ( |3.14| ) , where A may now be interpreted as the length of the world line. 
This restriction is usually assumed to mean that the length of the world line is the only 
remaining parameter after gauge-fixing. What we find instead is that there is a discrete 
infinity of classes, each with one or more continuous parameters. One member of this class 
comes from configurations A(r) which are everywhere positive; this class can be gauge-fixed 
to have 

A(r) =0, / drA(r) = A (3.15) 



with < A < oo. Similarily, the set of everywhere negative A(r) can be gauge-fixed as in 
( |3.15| ) but with — oo < A < 0. Keeping the first class alone gives the Feynman propagator 
for particles, while the second gives its complex conjugate.0 But we also have for example 
the class of A(r) which are positive for < r < n, negative for n < r < 1. The group of 
orientation preserving diffeomorphisms can gauge fix this to 

A(r) = for r ^ n, f rfrA(r) = Ai, f rfrA(r) = A 2 (3.16) 

JO Jt x 

with < Ai < oo, — oo < A2 < 0. We would like to identify this sector as the contribution 
to the amplitude to start with a state of type 1 and end with a state of type 2 (the 
off-diagonal element D12 of the matrix propagator). Similarily, all sectors beginning and 
ending with A > (thus having an even number of changes of the sign of A) contribute to 
D\\. We can add together these sectors for Dn once we choose the factor to be attached 
to each change in the sign of A. Choosing this factor is equivalent to choosing the regulator 
matrix M, and an explicit summation of sectors reproduces the matrix propagator. 



1 A restriction to the range (0, 00) for A can be naturally obtained using a Newton- Wigner 
formalism [13|. Here the particle travels only forwards in the time co-ordinate X°, thus it is not 
a co- variant approach. 
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3.3. BRST formalism. 

One might wonder if the more formal BRST quantisation of the relativistic particle 
would resolve the above issues about different possible quantisations. We follow the nota- 
tion in |jj4| . We introduce the canonical conjugate tv for A ([A, tv] = i) and ghosts (77*, Pi), 
(t) 2 ,V2) ([v l ^j] = f° r the two constraints tv = and \{P 2 — m 2 ) = respectively. 

The BRST charge 

2 

Q = ^tv + "L^-m 2 ) (3.17) 
is nilpotent (Q 2 = 0), and gives the variations: 



SX^ = -i[X», Q] = rfP^ SPf, = 

5X = r] 1 5tv = 

5?] 1 = 6V1 = -it ~ 

dr] 2 = 5V 2 = ~\{P 2 -m 2 ) 



(3.18) 



The equation of motion (obtained after gauge fixing) gives r] 2 = rj , which agrees with 

But we can define another nilpotent BRST charge 

Q> = ^'tv'X' + \' 2 {P 2 -m 2 ) (3.19) 

which generates the symmetry 

5X»' = rf'X'P^' 8PJ = 

5X' = X'r] 1 ' 5tv' = — 77 1 tx' ~ 

5r] x ' = 5V[ = -AV ~ 

Sri 2 ' = 8V' 2 = -\{P 2 ' -m 2 )~ b 

2 



(3.20) 



The equation of motion gives rj 2 ' = rj 1 A', which suggests that we should identify the above 
symmetry with S 2 - 

The symmetries Q and Q' are related through the identifications 



TV 



tt'A', A = logX' (3.21) 



all other primed variables equalling the unprimed ones. From ( |3.21|) we find that — 00 < 
A < 00 corresponds to < A' < 00. If we perform a path integral with the primed variables 
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and sum over both positive and negative A' then we are summing over more than is being 
summed in the unprimed variable path integral. 

We thus see sources of ambiguity on the quantisation of the relativistic particle work- 
ing with a Fadeev-Popov approach in sec 3.2 and a BRST approach in sec. 3.3. In fact 
the action we start with, ( |3T5| ) , is itself ambiguous because of the two possible signs of the 
square root. The particle trajectory would keep switching in general between timelike and 
spacelike, and at each switch we have to choose afresh the sign of the real or imaginary 
quantity obtained in these two cases respectively. This suggests that the world line con- 
figuration should be described by the pair {X M (r), cr(r)} with a = ±1 giving the choice of 
root. Evaluating the quadratic form of the action (given in ( |3.9| )) classically we find the 
sign of A to be related to the sign of the square root chosen for ( p.5|) . 

The above discussion suggests a close connection between the ambiguities found in 
three different approaches to the quantum relativistic particle, and it would be good to 
determine if they indeed are the same. For the rest of this paper we simply adopt as 
basic the picture of two complex conjugate propagations on the world line, with switching 
between them possible through the regulator matrix. 



4. A curved space example: spacetime with expansion. 

Consider the free scalar field ((2J.) with A = 0) propagating in 1 + 1 spacetime with 
metric 

ds 2 = C(r])[dr] 2 - dx% -oo < 17 < 00, < x < 2n (4.1) 

C{rj) = A + B tanh k??, A > B > (4.2) 

The conformal factor C{rf) tends to A ± B at r\ —>■ ±00. The limit k —>■ 00 gives a step 
function for C(rj); the Universe jumps from scale factor A — B to A + B at 77 = 0. We will 
work in this limit to ensure simpler expressions. 
For r\ — f — 00 it is natural to expand (p as 

00 

<j>(ri,x) = V -= ; (a n e-"^ +mx + a\e lw ^-™) (4.3) 

^ • / 2n y /2uj^ 



n — =00 



with 

uj~ = (n 



2 + (A- B)m 2 ) 1/2 >0 (4.4) 
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We define the 'in' vacuum by 



&n|0 >in= 0, for all n (4-5) 

Similarily, for 77 — > 00 we write 

00 

0(t/, x) = V = ; (a n e- i < r ' +inx + 4e <w »"- in!B ) (4.6) 

with 

w + = (n 2 + (A + B)m 2 ) 1/2 >0 (4.7) 
The 'out' vacuum is defined through 

a n \0 >out= 0, for all n (4-8) 

The 'out' vacuum does not equal the 'in' vacuum, even in the free theory: 

n>0 



|0 > out = Coe^ a ° a ° n C n e*»< a -«\0 > m (4.9) 



Wit + UJj 



C n = (l-6 2 ) 1 / 2 , n>0. (4.10) 



The Bogoliubov transformation (|4.9| ) says that an 'out' observer will find particles in his 



frame as r\ — * 00, if the 'in' observer sees a vacuum. Since the Bogoliubov transformation is 
given by the exponential of a quadratic in the field, we guess that our formalism developed 
in the preceeding section should apply. In other words, the effect of the flux created by 
spacetime expansion can be incorporated by a change in the regulator matrix of the first 
quantised path integral. We demonstrate this explicitly in our example. 

The first quantised path integral gives @ (we denote the pair (77, x) by z) 

<z 2 \ I™ d\D[X]e-^>^ x2 ' X+ ™ 2x) \ Zl > = ^<0TOW(^i)]|0>in (4n) 

JO out < > in 

(We have integrated out p in the phase space path integral.) It is not surprising that both 
vacuua appear in this quantity; after all the action and measure are covariantly given and 
should not distinguish the past or future as special. Using — i in place of i in the exponential 
gives in < 0\f[4>(z 2 )(p{zi)]\0 > out /in < 0|0 > out . (T denotes anti-time-ordering.) 

In studying kinetic theory we typically wish to specify the density matrix in terms of 
the 'in' states, e.g. p = Y1 e~^ n |n >i n j n < n\/(^2 e~^ n ) where \n >; n gives the occupation 
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number n state for some positive frequency mode at past infinity. For (3 — > oo, p = p = 
|0 >in in < 0|. We use po for our illustration; it should be straightforward to consider 
both expansion of spacetime and an initial exponential distribution of particles, by putting 
together the results of this section and the last section. 

To develop a perturbation theory using po we need a real time contour running from 
rj = — oo to 77 = oo, and then back to rj = — oo where we insert po and take a trace to close 
the path. This perturbation theory requires a matrix propagator D: 



n(y r x / in < O|T[0(z 2 )0(zi)]|O >in in < 0\^(z 1 )^)(z 2 ) |0 > in \ . . 

i; ' 2,; ! ' 1 in < O|^ 2 )0(^)|O > in to <O|f[0(^)0(«i)]|O> < 1 {Ll - ) 



Using the 'matrix action' i(— U — m 2 )o~3 — el in the path integral ( [4.1 1|) gives the matrix 
propagator D° = diag{ out < 0\T[(fi(z 2 )(j)(zi)} |0 > in , in < 0\f[(fi(z 2 )(j)(zi)}\0 > out ] We wish 
to change only the operator multiplying e, and we wish to get (|4.12| ). 

Let us set up the calculation of Green's functions in the first quantised formalism. We 
need eigenfunctions of the world line Hamiltonian: 

(□ + m 2 )i/} s (r),x) = -Hi(; s (r},x) = si/j s (rj,x) (4.13) 

The following is a complete set: (— oo < n < oo) 

,2 „2 



,2 



n o n 



m + — — < s < m + 



A+B A-B 

^inx 



'2-k 



= [-(! + -^)e-™^ + -(l--^Y v -% V <0 



V27r 2 2 
z> + = [-(A + B)(m 2 -s) -n 2 } 1/2 > 0, v_ = [(A - B)(m 2 - s) + n 2 ] 1/2 > (4.14) 



i n2 



— oo < s < m + 



A + B 



V Z7T 



r "" rl (l + ^e" 1 ^ + i(l - ^)e lv ~% rj < 



V27T 2 2 
(A ± B){m 2 — s) + n 2 , sign ( v + )— sign [v-] (4.15) 
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These functions are normalised as 



(ipi>'n',^v + ,n) = / dr)dxC{rj)\l)l, n ,(r],x)ip i)+ ^(r],x) 



= 5 nl>n Tr{A + B 



v 2 _ + v 2 



(4.16) 



2z/_z>_|_ 



( \ 2 2 2 

= V.n^ + ^H o + 8{i/ + -u+) + U -~ U+ d(u' + + u+)} (4.17) 
Let us first recover the propagator (|4.11] ) in this formalism. Let rj',r) > 0. The range 



2 2 

m 2 + iTo < s < m 2 + 



A—B 



gives the contribution 



V/2B/(A-B)| 



E 



f+ ,Tl I 



< Vtf, nk&> 5(r7+ -v+)[{A + B) 



TV 



z/ 2 + P 2 , 



4- 2tt y * 



2z/_z/_|_ 



^_ + v\ iy\ + n 2 + m 2 (A + S)) 

(4.18) 



o 2 

Similarily the range — oo < s < m + J^_ B provides the contribution 



E 



i(x' — x) 

~2^ 



— \e iu+( - v ~ v) + 

+ 27T L U + + V 



(z/ 2 - n 2 - m 2 (A + S) + ie) 



(4.19) 

There is a branch cut in the complex v+ plane joining 1/4. = ±zy^P^|n|. The 1/4. integral 
in (|4.19|) has a discontinuous jump across this cut for the part multiplying e lv+ ^ 11 +r] \ The 



contribution from (|4.18|) can be added to (|4.19|) , however, with the identification z?_|_ = iv+. 
This results in a contour passing over the cut. Evaluating the resulting contour integrals 
one obtains the result 

gin(x'— x) -y 



2tt 2ut ' 



ibJ+\v'-v\ 



+ 



1 w~ n 

n 



e -iw+(u+»j) j for v r] > > o 



(4.20) 

which may be readily verified in the operator language using ( |4.6| ) and ( |4.9| ). 

In the flat space examples of section 3 energy-momentum conservation implied that 
the regulator matrix was diagonal in the basis of eigenf unctions of — (Cl + m 2 ) given by fixed 
(po,p). In the time-dependent situation that we have now, energy is not conserved. The 
regulator matrix acts within the eigenspace of a fixed value of s. Thus for each n we need 
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to consider a 4x4 matrix, which acts on the column vector ({f} n , fs, n } + i {fs,m fs,n\~)- 
Here the first pair of functions propagate on the world line as e~ ls while the second pair 
propagates as e ls . Within each type (+ or — ) we have two linearly independent functions 
for any given s near the mass shell s = O.i Since we wish to compute the propagator with 
the 'in' vacuum density matrix, we choose in this space the basis which at the mass shell 
becomes {ui~ > 0) 



f 1 

J n 



T) < 



fn(v,x) = f n *(r},x) 



V27T ^ UJn 
1*, 



77 + -(1 



rj > 



(4.21) 



We want the matrix propagator for the density matrix p = |0 >i n i n < 0| and the 
time path starting at t = — oo, going to t = oo and then back to t = — oo. A calculation 
similar to the above yields that the propagator (|4 . 1 2[) is obtained if the regulator matrix 
in the space of the nth Fourier mode e Lnx is 



M 





1 














1 




-2 






1+4B 




-2 




1 






1 


{ 





1+ S B;i 



/ 



(4.22) 



where 



B n = — — 



2 Un + OJr. 



(In this computation we need to note that lim e ^o~ 



7r5(x), but lim 



(4.23) 
= 0.) 



5. Discussion. 

We have taken the view that in a theory of quantised matter and gravity the back- 
ground for perturbation theory should not only be a specification of classical values of 
fields, but also a specification of 'exponential of quadratic' particle fluxes. We know that 
these two different aspects of the background arise naturally in flat space kinetic theory. 
With curved space, it becomes natural to consider the kinetic theory and to not construct 

2 The regulator matrix needs to be defined only in the infinitesimal neighbourhood of the 
mass shell. 
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a 'vacuum' theory at all. The reason is that starting with physically acceptable conditions 
in the past, say, particle fluxes can be created in the co-ordinates natural in the future. 

The usual first quantised approach of string theory, applied to the scalar particle, 
gives a 'vacuum' theory, where certain 'in-out' vacuum boundary conditions are chosen at 
temporal infinity. These boundary conditions appear unnatural for physical purposes, so 
we would like to be able to move to a more general class of states at the boundary. In 
particular we would like to be allowed a radiation flux at the past time boundary, as in 
the radiation dominated Cosmologies. We find that in the first quantised language there is 
a natural way to obtain this more general theory. Quantisation of the relativistic particle 
indicates that we should consider both propagations e~ lHX and e lHX on the world line. 
The regulator matrix needed to regulate this path integral need not be diagonal in these 
two modes of propagation. Off-diagonal terms encode an 'exponential of quadratic' density 
matrix and a choice of time path for perturbation theory. 

One special case of particle flux, the flux for constant temperature may be studied 
without the real time formalism. One studies the theory on spacetime with time rotated to 
Euclidean signature and identifies t with t—i/3 [|15|. But constant temperature is unnatural 
in a theory with gravity, as the particle density gives a gravitational field, which gives an 
acceleration of the scale factor of the Universe (unless we carefully balance the matter 
density with a Cosmological constant). The changing scale factor violates the constant 
temperature requirement needed for the Euclidean time trick to work.i 



Leblanc |16| studied the real time formalism for open and closed strings, for the case 
of constant temperature. The propagator was computed in the 'thermo-field dynamics 
language, which used the Niemi-Semenoff time path, and so was given by (|3.2|) , (|3.3| ) . For 
this time-independent situation amplitudes were computed and the Hagedorn temperature 
recovered. In [DJ The imaginary time formulation was used to compute the Background 
field equations for the closed string at constant temperature. 

We should distinguish two different limits in which the physics of fluxes may be stud- 
ied. One limit is where the collisions are so rapid that approximate thermal equilibrium 
is maintained at all times, and we need only let (3 be a function of time. The other limit 
is that of kinetic theory, where we assume that collisions are rare; particle wavefunctions 



3 More general flux situations may be obtained by making a canonical transformation on the 
variables living on the world line and compactifying the new 'time' co-ordinate after analytic 
continuation to imaginary values. But such a description does not appear physically illuminating. 
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evolve on the time-dependent background, and collisions between these particles are taken 
into account by perturbation theory. Our approach assumes the latter limit. 

The limit e — > implies that the effect of the regulator matrix M is felt only on-shell 
(i.e. for p 2 — m 2 = 0). Equivalently, we may say that only world lines of infinite length 
(A = oo) see the regulator matrix. To see this, let M and M' be two different regulator 
matrices. We write 

L _ poo 

HX-eMX , / j\ -iHX-eMXi 



pLi _ pOO 

D M > = lim£_ K »lim € _> 0+ {[ / dXe~ iHX - eMX + / dXe 

Jo Jl 



poc _ _ poo 

+ [ / d\e- iHX -* M ' x - / d~Xe- lHX -" MX ] (5.1) 
Jl Jl 

+ [[ d~\e- iHX - eM ' x - [ L d\e- lHX -" MX ]} 
Jo Jo 

The first square bracket on the RHS is Dm, the last vanishes with the indicated limits, 
while the second has support only on world lines of infinite length. 

The analogue of the above statement for strings is that the effect of p is to give a 
contribution to the boundary of the moduli space of Riemann surfaces, where a homolog- 
ically trivial or non-trivial cycle is pinched. (It is important to have Minkowski signature 
target space, and correspondingly a Minkowski signature world sheet, to allow the on-shell 
condition for the particle flux.) A /^-function calculation for the string world sheet theory 
would have to take into account such pinches while considering the small handle contribu- 
tion studied by Fishier and Susskind |18[] . This calculation should yield a relation between 
the classical fields and the particle fluxes, rather than just among the classical fields giving 
the background. 
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Appendix A. Wick theorem for 'exponential of quadratic' density matrices. 

We wish to establish Wick's theorem for density matrices of the form 

p = e aaV e -0ato e 7oa ( Al ) 

A string of creation and annihilation operators can be brought to normal ordered form 
in the same way as for the usual Wick theorem in the vacuum. What we need to show in 
addition is that 

ifr{pat...ata...a} = r% tflt },,,%M J^W (A2) 
Trp lP J ^ Trp Trp Trp y ' 

where the RHS has a summation over all possible pairings of the a^, a operators on the 
LHS. We sketch below some of the steps involved in the derivation. 
Note 

gaa 1 a f e ~/3a 1 a _ g - (3a 1 a g a' V a 1 " (A3) 

with a' = ae 2 ' 3 . A straightforward calculation gives 



4«7 

j.rp = <± - e ' ) [i - — 

which we may rewrite as 



^ = (i-^rii-j^^r 1 ' 2 (a.4) 



Trp = TVe-/3 ata e a ' ata V aa } = (i_ e -£)-i[i- 4a \xJ ~ 1/2 (A.5) 

(1 — e ^) 2 

From ( |A.1|) we see that in (|A.2| ) there must be either an even number (2p) of 'a' oscillators 



and an even number (2q) of 'a' ' oscillators, or an odd number (2p + 1) of 'a' and an odd 
number (2q + 1) of 'a^' oscillators. Assume first that we have the former case. Then the 
LHS of (|A.2|) is obtained as 

Tv{p(a^\a) 2p } = -L(d a ,)«(d^ p Trp (A.6) 

Trp 

(Partial derivatives are taken with a', 0, 7 as independent variables, unless otherwise 
mentioned.) In particular, 

< aV > = -^da'Trp = 2j e - 2(3 /K = A 
Trp 

<aa> = — <9 7 Trp = 2a'e- 2f3 /K = B (A.7) 
<a+a> = - J-df}[Trp] a „ = e-P(l-e-P)/K = C 
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where in computing C, dp is a partial derivative with a, 7 held fixed. Here 



K = (1 - e~ p f - 4c/ 7 e" 2/3 (A.8) 

We find 

Trp = eP[C 2 -AB] 1 ' 2 = K~ x l 2 (A.9) 

For fixed (3 



d a <A 


= 2A 2 


d 1 A = 


2c 2 


d a ,B 


= 2C 2 


a 7 s = 


2s 2 




= 2AC 


a 7 c = 


2BC 



(A.10) 

Using the above formulae, we can establish (|A.2|) by induction. Suppose (|A.2|) holds 



with 2p operators 'a' and 2q operators 'a^'. A typical term on the RHS would have the 
form FA ni B n2 C n3 , where F is a constant and n±, 77,2, ns > 0. To establish the result for 
2p operators 'a' and 2q + 2 operators 'a'' we get for the LHS of ( |A.2| ): 

d' le^lC 2 - AB) 1/2 FA ni B n2 C ns ] = F\A ni+1 B n2 C n3 

eP(C 2 -AByP al ■ { } [ ( A-11 ) 

+ 2n 1 A ni+1 B n2 C n3 + 2n 2 A ni B n2 - 1 C n3 + 2n 3 A ni+1 B n2 C n:i } 



The first term on the RHS of ([A. 11 ) gives the pairing of the two new operators with each 



other. The second term gives the n\ ways to choose an existing pair (aW) and to contract 
the new operators with members of this pair instead. The third term corresponds to 
choosing an (00) pair in the original expression and contracting the 'a' operators with the 
new 'aT' operators instead. The last term corresponds to exchanging the in an existing 
a^a pair with one of the new operators. It is easily seen that this generates all the new 
terms required on the RHS of ( |A.2| ) for the induction to hold. 

To work with the case of an odd number of a and operators we start with the 
expression CTvp + 2^d 1 Trp = Tr{e~^ ata e a ai ai ae iaa } in place of Trp, and proceed as 
above to introduce extra aW and aa pairs in the induction. 
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